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Abstract 


We present a new analysis of observations of comet 29P/Schwassmann-Wachmann 1 taken on UT 2003 November 

21, 23, and 24, similar to a previous investigation of the observations (Stansberry et ah 2004), but using the most recent 
Spitzer data pipeline products and intensive image processing techniques. Analysis of images from the IRAC 5.8 & 8.0 
/im bands and the MIPS 24.0 & 70.0 pm bands resulted in photometry measurements of the nucleus after a suite of coma 
modeling and removal processes were implemented. SWl was not identified in the 5.8 pm image from the previous work so 
its incorporation into this analysis is entirely new. Using the Near Earth Asteroid Thermal Model (Harris 1998) resulted 
in a nucleus radius measurement of R = 30.2 km and an infrared beaming parameter value of 77 = 0.99 lo!i9* 
measured an infrared geometric albedo, 795.8 = 0.5 ± 0.5. Extrapolating a 0.04 V-band albedo and using a normalized 
reflectivity gradient S' = 14.94 ± 1.09 [% (1000 A)“^] (Duffard et al. 2014) we recover an infrared albedo of 795.8 = 0.31 
in the near infrared consistent with the value recovered from thermal modeling. The dust composition extracted from 
IRS spectra are very comet-like, containing mainly amorphous ferromagnesian silicates (but with a minority of crystalline 
silicates as well), water ice, and metal sulfides. 
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1. Introduction 


Comet 29P/Schwassmann-Wachmann 1 (SWl) is a unique comet and excellent observational target to investigate for 
improvements in our understanding of cometary nuclei and their activity drivers. SWl is in a nearly circular orbit with 
eccentricity e = 0.041830, semi-major axis a = 6.009674 AU, and inclination i = 9.3763° (orbital elements from lAU 
Minor Planet Center, Epoch 2014 Dec. 9.0 TT = JDT 2457000.5). The orbital properties of SWl also classify it as a 
Centaur (Jewitt 2009). With a nearly constant background level coma (believed to be driven by the supervolatile CO, with 
observations of CO productions rates made by Eestou et al. (2001), Gunnarsson et al. (2001), and Paganini et al. (2013)), 
SWl undergoes frequent outbursts in activity with an increase in magnitude ranging over several orders of magnitude 
(Gronkowski 2014; Kossacki and Szutowicz 2013; Trigo-Rodriguez et al. 2010). The combination of SWl’s orbit and 
nearly constant insolation raises the question of what activity drivers are causing the outbursts. With a seemingly stable 
thermal environment why would the nucleus experience frequent excursions from a steady state of activity? This puzzling 
behavior is why SWl is an excellent target for continued study and analysis. 

The radius distribution of the Gentaur population is still an area of active research. Due to their large heliocentric 
distances, observational size measurements are challenging. A radius range of 2 km to 41 km has been observed for 
Gentaurs by Jewitt (2009) and the first discovered Gentaur 2060 Ghiron has a 210l}J km radius (Eornasier et al. 2013). 
Where does SWl he in this distribution? Several groups past and present have studied SWl in order to measure properties 
of the comet’s nucleus. Due to the continuous activity of SWl producing a persistent coma, direct measurements of the 
nucleus are currently unattainable. Eurther, observational measurements of SWl are hindered by its large heliocentric 
distance. Radius measurements of SWl show a large range: 20.0 ± 3.0 km (Cruikshank & Brown 1983), 8.6 ±0.1 km 
(Meech et al. 1993), 27.0 ± 5.0 km (Stansberry et al. 2004), km (Stansberry et al. 2008). This large spread of 

radius values shows our lack of knowledge of basic physical properties of SWl’s nucleus. Using Spitzer thermal images 
we have made an attempt at constraining the nucleus’ effective radius, infrared beaming parameter, and albedo. Using 
Spitzer spectral observations, we have also constrained the coma’s dust composition. 
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Table 1 


SWl Observation Summary 


Instrument 

AORKEY 

Epoch^ 

h 

Band^ 

t ^ 

^exposure 


IRAC 

6068736 

21, 06:03 

5.8 

150 

5.51 




8.0 

150 


IRS 

6068992 

23, 07:18 

SLl 

122 

5.54 




LLl, LL2 

59 


MIPS 

7864064 

24, 15:05 

24.0 

66 

5.56 




70.0 

40 



^ Day of 2003 November, UT at beginning of observations. 

^ Imaging wavelength (fim) or spectral band. 

^ Total exposure time per pixel for each image frame, seconds. 

^ SWl-Spitzer distance (AU). Heliocentric distance was 5.73 AU, solar phase angle was 10.0° and the tracking rate was 8".6 hr“^. 


The organization of the paper is as follows. Section 2 gives an overview of the Spitzer observations, including details 
pertaining to improvements of the Spitzer BCD pipeline and emphasizing the improvements of the presented results over 
their 2004 counterparts. The coma modeling and removal procedure is described in Section 3 along with an application of a 
Near Earth Asteroid Thermal Model (NEATM, Harris 1998) to the extracted nuclear infrared photometry measurements. 
In Section 4, results from IRS spectrum analysis are presented. Section 5 gives an overview of the results and their 
implications on the thermal evolution experienced during the dynamical evolution of SWl. 


2. Observations 

Spitzer observations of SWl were acquired on UT 2003 November 21, 23, and 24 during the In-Orbit Checkout 
(IOC)/Science Verification (SV) period (Werner et al. 2004). Table 1 gives the details of each observation (similar to 
Table 2 from Stansberry et al. 2004). All three instruments on Spitzer were used for the observations: Infrared Array 
Camera (IRAC) (Eazio et al. 2004), Infrared Spectrograph (IRS) (Houck et al. 2004), and Multiband Imaging Photometer 
for Spitzer (MIPS) (Rieke et al. 2004). The details of each observation are given in the following three subsections. 


2.1. IRAC 

The Infrared Array Camera (IRAC) is a multi-channel infrared camera capable of imaging in 3.6, 4.5, 5.8 and 8.0 
/im bands. A detailed description of the IRAC can be found in the IRAC Handbook (IRAC Instrument and Instrument 
Support Teams 2013). Each band images a field of view (EOV) of 5.2'x 5.2' and a pixel scale of 1.2"/pixel. 

The IRAC was used in high dynamic range (HDR) mode during the observations. Unfortunately, a star was close to 
SWl in the EOV producing a diffraction spike in the images. This rendered the 3.6 /im and 4.5 /im observations useless 
other than to acquire upper limit photometry. Of importance is the identification of SWl in the 5.8 /im band, which had 
not been used in the previous analysis (Stansberry et al. 2004). The 5.8 and 8.0 /im observations entailed acquiring five 
30 second exposures of SWl. The software MOPEX (Mosaicking and Point-source Extraction, Makovoz et al. 2012) was 
used to calibrate and stack the basic calibrated data images (BCDs) acquired from the Spitzer Rentdige Archive, generated 
from version S18.25.0 of the Spitzer pipeline. The final stacked, calibrated, and cropped 5.8 /im and 8.0 /im images are 
shown in Eigure 1. Each image has an effective exposure time of 150 seconds. 

In the 5.8 /im image, SWl was identified close to the mentioned diffraction spike. A broad diffuse coma was not 
observed in this band. On the other hand the 8.0 /im image shows traces of a diffuse coma in the northeast and southeast 
region of the image. The diffraction spike is also in this region, but was suppressed using correction techniques described 
in the IRAC Handbook. 
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(a) (b) 


Figure 1: (a) 5.8 /iui and (b) 8.0 fim combined and calibrated images centered on SWl. Each has a FOV of 1.01' x 1.01' 
and equatorial north is up and east is to the left. The solar direction is identified by the yellow arrow and SWl’s projected 
heliocentric velocity by the red arrow. The scale bar at the bottom right of (a) shows a distance of 20" and the circle in 
the top left of the images represents the size of the PSF’s first minimum. A star’s diffraction spike artifact is clearly seen 
in the 5.8 /im image and partially visible in the 8.0 jam image. 


IRS 

Spectra of SWl covering 7.4-38.0 /im were obtained with the IRS operating in the Long-Low order 1, Long-Low order 
2, and Short-Low order 1: LLl, 19.5-38.0 /im, LL2, 14.0-21.3 /im, and SLl, 7.4-14.5 /im. The angular size of the slits on 
the sky are 11" x 170" for the LL modes and 3".6 x 57" for the SL mode. The exposure time for each of the observations 
can be seen in Table 1. A detailed summary of IRS and spectral reduction and analysis techniques can be found in the 
IRS Instrument Handbook (IRS instrument Team & Science User Support Team 2011). The Standard Staring Mode of 
operation was used, which included a nodding of the telescope. 

2.3. MIPS 

The Multiband Imaging Photometer for Spitzer (MIPS) is an infrared imager with 24.0, 70.0, and 160.0 /im bands. A 
detailed summary of MIPS can be found in the MIPS handbook available through the Spitzer Heritage Archive (MIPS 
Instrument and Instrument Teams 2011). SWl was observed in all three bands, but only recovered in the 24.0 /im and 
70.0 jam images. Observations with MIPS involved scanning a large region around SWl, which resulted in most pixels 
in the mosaic having an ^70 second exposure for the 24.0 jam and ^40 second exposure for the 70.0 jam image. These 
images, after MOPEX mosaicking and calibration, can be seen in Figure 2, each having an 8.0'x8.0' FOV. Figure 3(a) 
shows the full 24 jam mosaic of SWl along with intensity contours showing the coma shape and structure. This figure can 
be related to Stansberry et al. (2004) Figure 1, but we have oriented the images with Equatorial North up as opposed to 
Ecliptic North. Also, Figure 3(b) shows the 24 jam image with a l/p profile removed similar to Stansberry et al. (2004) 
Figure 2 lower left panel {p is the skyplane projected cometocentric distance). 

MIPS mosaicked images are notorious for image artifacts in the final MOPEX generated mosaic. Each of these image 
artifacts are described in detail in the MIPS handbook, along with correction algorithms to mitigate their impact on the 
final mosaic. Processing for the 24 /am image involved applying a median time filter per pixel, which removed the long 
lived dark latent artifacts caused by a > 50 Jy source being imaged after the last array annealing, but before the SWl 
observations. The 70 pm image suffered from artifacts resulting from slow response variations on the array leading to 
streaks and stim latents due to the stim flashes used to measure the slow response of the array. Two correction algorithms 
were implemented on the BCDs before mosaicking which improved the quality of the final image: high-pass median time 
filter per pixel and column median value subtraction. Traces of a “jail-bar” artifact can still be seen in Figure 2 (b). 
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Figure 2: (a) 24.0 fim and (b) 70.0 fim mosaicked and calibrated image centered on SWl. Each has a FOV of 8.0' x 8.0' 
and equatorial north is up and east is to the left. The solar direction is identified by the yellow arrow and SWl’s projected 
heliocentric velocity by the red arrow. The scale bar at the bottom right of (a) shows a distance of 100" and the circle in 
the top left of the images represents the size of the PSF’s first minimum. 


3. Nuclear Properties 

Measuring nuclear properties of active comets is hindered by the presence of emission from the gas and dust contained 
in the coma. For an accurate assessment of the thermal radiation emitted by SWl’s nucleus the coma flux must be 
modeled and removed. This process is one of painstakingly and meticulously comparing the residuals of an image, with 
a model coma removed, to a point-spread function (PSF) of the optical system used for the observations. It should be 
pointed out that for most comet observations the nucleus is unresolvable in the image and its contribution to the image 
is a PSF. Thus, using a scaled PSF as a measure for the flux contribution is warranted. Our software was developed by 
our team and is based on the technique described by Lamy and Toth (1995) and used many times since then (e.g., Lamy 
et al. 2006, 2011; Fernandez et al. 1999, 2013; Kelley et al. 2013), in particular, to correctly predict the sizes of nuclei 
before their flyby by spacecraft (Lamy et al. 1998; Fernandez et al. 2003; Lisse et al. 2009). In this section a description 
of the nucleus extractions techniques is discussed, the resulting nuclear spectral flux density measurements are presented, 
and nuclear properties extracted from thermal modeling are shown. 

3.1. Coma Removal 

The modeling and removal of the coma from a comet’s image is not a straightforward procedure or “one size fits all” 
technique that can be applied to comet observations. Each image must be analyzed to determine a modeling and removal 
method best suited for a successful point source extraction. Two modeling techniques were found to be necessary for the 
Spitzer images: a scaling of a system PSF for the 5.8 jam and 8.0 jam images and a more complex coma modeling routine 
for the 24.0 jam and 70.0 jam bands. Each method will be described in detail in the following paragraphs. Both have been 
tested with synthetic comet images to verify coma removal capabilities. 

For the 5.8 fim and 8.0 /am images, shown in Figure 1, there was not sufficient flux from the coma to produce a usable 
coma model. In addition, the removal of the diffraction spike from the images resulted in an artifact that interfered with 
any coma flux present. This method involved first finding the pixel location of the comet’s centroid and choosing this to 
be the location of the nucleus. Next, a PSF was centered on this location and scaled to match the level of the comet’s 
flux. Initially, a PSF generated from the program STINYTIM (Krist 2006) was used for this procedure, but it was found 
that this PSF did not represent well the structure observed in the comet’s image or field stars also in the image. One of 
the field stars from each band, sufficiently distant from the diffraction spike, was used to represent the PSF. SWl was 
moving slow enough during each image’s exposure that the projected distance traveled was less that a pixel scale length, 
which was approximately 4800 km/pixel. Once the PSF was scaled it was subtracted from the comet image, resulting 
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Figure 3: (a) 24 fim MIPS mosaic of SWl showing contours to highlight the coma and tail shape. The mosaicked image 
has dimensions of (17' x 52'). The contour levels of the image are (0.2, 0.5, 1.0, 2.5, 5.0 and 10.0) MJy/sr. (b) 1/p removed 
24.0 yum image similar to Stansberry et al. (2004) Figure 2 {lower left) panel, which shows enhancement in the south-east 
region produced by the single jet observed in the 2004 analysis. It should be noted that the jet feature is due to solar 
radiation pressure and not from rotation as was previously believed. The scale and orientation of this image are the same 
as Figure 2(a). 
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in what should be only flux from the coma and background. A best-fit PSF was taken to be one that minimized the 
standard deviation of the residual of pixels from a region centered on the comet’s centroid. For the 5.8 /im image there 
was not much come and the diffraction spike also presented a problem so this was the best approach found to measure 
the nuclear flux. This scaled PSF was then taken to represent the nucleus’s contribution to the image flux and used 
for photometry measurements. Figure 4 shows radial profiles of the 5.8 /im comet image and scaled PSF for a selection 
of azimuthal angles. Figure 5 shows the 5.8 /im comet image, scaled PSF image, and residuals after PSF subtraction. 
Similarly, Figures 6 and 7 are for the 8.0 jam image. 
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Figure 4: Radial cross-section of the 5.8 /im image showing the scaled PSF. The angular descriptor indicates the position 
angle of the radial cross-section. The coma and artifact contributions to the images can be seen in the comet profiles. 
Note that there is not much coma in the core. 
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Figure 6: Radial cross-section of the 8.0 /iui image showing the scaled PSF. The angular descriptor indicates the position 
angle of the radial cross-section. Note there is not much come coma present in the core. 



(c) 


Figure 7: The (a) 8.0 jam comet image, (b) the scaled PSF, and (c) the residual after PSF subtraction. Notice in (c) 
there is an over subtraction to the south of the centroid and under subtraction to the north of the centroid. This can be 
attributed to asymmetry in the star generated PSF and SWl nucleus image. 


The 24.0 jam and 70.0 /im images had sufficient coma present for a more elaborate coma modeling routine to be 
implemented on these images (Lamy et al. 2004; Lisse et al. 2009; Fernandez 1999). The procedure for this method is to 
first take azimuthal pie sections of the comet image, centered on the comet’s centroid. Radial profiles of each pie section 
are generated. The angular width of the sections was a variable for the procedure and was selected independently for 
each image. The S/N of the image determined the minimum angular spread usable, with smaller angles being desired for 
higher fidelity in the model coma. Again, since the nucleus is unresolvable to Spitzer at this distance it will only contribute 
fiux in the region of the PSF. In theory the nucleus should only image as a PSF on the image plane. Thus, any fiux 
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Figure 8: Shown are radial profiles of the 24.0 /im comet image, STINYTIM PSF, and synthetic coma model. The PSF 
effectively goes to zero at a certain radial position which can be seen to be around 30 pixels. Comet flux beyond this 
region was used for the fitting procedure. Extrapolation of the coma model to the center shows the excess contribution 
to the comet image attributed to the nucleus. 


present outside of the region of the PSF should be contributions from the coma and background. An example of these 
radial profiles from both the 24.0 /im and 70.0 /im image and their corresponding PSFs can be seen in Figures 8 and 9. 
If regions radially distant from the nucleus (regions outside the dominance of the PSF) are chosen, a coma model of the 
form can be fit to the profiles, generating a synthetic coma. The parameter 0 is the position angle (PA) of the 

pie section used for the radial profile. This procedure was implemented on each of the azimuthal pie section producing 
the coma model shown in Figure 10(b) for the 24.0 jam band and Figure 11(b) for the 70.0 /im band. The coma model is 
then subtracted from the comet image yielding the PSF’s contribution to the image. A comparison of the residual PSF 
and an STINYTIM-generated PSF is shown in Figures 10,11 (c) and (d). Close inspection of the two shows a high degree 
of similarity (the first order bright fringe as an example), giving a level of validation to our modeling technique. 

3.2. Nuclear Photometry 

Once the coma was successfully removed from the images, photometry was used to measure the spectral flux density in 
each of the four observational bands. Since the coma removal process resulted in a scaled PSF representing the contribution 
to the spectral flux density from the nucleus, no aperture correction was needed for the photometry. Table 2 shows the 
photometry results after the application of color corrections. Uncertainties in the spectral flux density measurements 
are derived from analysis of the distribution of the flux measurements from the many coma models for each band. The 
procedure to arrive at the best-fit coma model involved varying parameter space of the model (centroid pixel location, 
angular size of the azimuthal pie section, size in pixels of the radial region used for fitting, center position of the radial bin 
of pixels) and observing how this minimized the difference between the residuals after coma removal for the case of the 
24.0 /im and 70.0 /im and the standard deviation of the central pixels for the 5.8 /im and 8.0 /im images. Each of these 
coma models was stored and the resulting residuals from each model was measured for the nucleus’ contribution to the 
spectral flux density of the image. The distribution of the flux measurements from the many coma models that we tried 
for each band was used to derive uncertainties for the photometry measurements. We varied the parameter space of the 
model (centroid pixel location, angular size of the azimuthal pie section, and radial region used for coma slope fitting) so 
as to determine how a model’s answer for the nucleus photometry depended on these input parameter assumptions. The 
scatter in the nucleus photometry was used to derive the uncertainties we list for the nucleus in Table 2. 

With the nuclear photometry measured, the spectral flux density measurements for a 9".0 radius aperture of the 
coma’s contributions were found for each of the IRAC and MIPS bands. These values can also be seen in Table 2. These 
values were found by subtracting the nuclear flux values from each image and then performing aperture photometry on 
the residual coma flux. Coma photometry is used to help tie together the orders that make up the IRS spectrum, and will 
be explained further in Section 4. The overall comet (coma + nucleus) photometry is well constrained, and the error bar 
for the coma is driven almost entirely by the error bar on the nucleus measurement. 
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Figure 9: Shown are radial profiles of the 70.0 /im comet image, STINYTIM PSF, and synthetic coma model, similar to 
the 24.0 jam profiles in Figure 8. 




(b) 
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(c) 


(d) 


Figure 10: (a) 24 /am comet image, (b) the synthetic coma model, (c) the residual after model coma removal, and (d) the 
STINYTIM 24 /am PSF. Notice the asymmetry seen in the comet image is modeled well in the coma. Also, the residual 
nucleus contribution has many features similar to the optical system PSF, indicating a good coma removal. 
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3.3. Thermal modeling 

The Near Earth Asteroid Thermal Model (NEATM, Harris 1998) is an improvement of the Standard Thermal Model 
(STM) in its applicability to many small bodies since it treats phase darkening in probably a more realistic way and since 
it lets the beaming parameter float as a free parameter. This model has been applied to many minor bodies (Mainzer et al. 
2011; Delbo et al. 2011) and was applied to the Spitzer photometry measurements, returning values for the effective radius 
(i?), IR beaming parameter (77), and geometric albedo (ps.g)- The IR beaming parameter is used as a fitting parameter 
in the NEATM, with a value of 77 = 1 applying for a perfectly memory-less, spherical surface with modest topography. 
Values of 77 > 1.0 imply lower surface temperatures than what would be observed for the ideal situation possibly caused by 
thermal communications between the day and night sides of the object. Topography and surface roughness can be inferred 
from values of 77 < 1.0. This is due to the observed surface temperature being higher than what would be observed from 
strictly applying the STM. Excess thermal flux can be thought to arise from emission from bowl shaped craters, where 
emission from the sides of the crater is absorbed at the bottom of the crater. This additional flux results in an increase in 
the equilibrium surface temperature when compared to the STM temperature. The beaming parameter is a way to adjust 
the equilibrium surface temperature distribution of the object and have it match the color temperature of the object’s 
observations. 

The 5.8 jam photometry value is too high to be explained by thermal emission alone. Inclusion of a reflected light 
component, i.e. a scaled solar spectrum, into our model however explained the 5.8 jam flux. This is what would have let us 
constrain the V-band albedo p, in principle, although it turned out that the uncertainties were too great for a meaningful 
result for p. Our best fit model (using minimization), making use of both thermal and reflected components, yielded R 
= 30.2 ^3 7 km, 77 = 0.99 ^0 26 ’ ~ ^ Figure 12 shows the measured spectral flux density values, NEATM 

+ Reflected curves, and best-fit spectral flux density values. Other parameters necessary for the thermal modeling are: 
bolometric bond albedo A = 0.012 (assuming a visible-wavelenth geometrical albedo 79 = 0.04 and phase integral relation 
q = 0.290 + 0.684G, Harris and Lagerros 2002 ), emissivity e = 0.95, and slope parameter G = 0.05 (the same assumption 
made by Eernandes et al. 2013 in their Spitzer survey of cometary nuclei). 
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Table 2 


Nuclear Photometry Measurements 


Band (//m) 

Coma Spectral Flux Density^ (mjy) 

9".0 Radius Aperture 

Nucleus Spectral Flux Density^ (mJy) 

5.6 

8.0 

1.8 ± 1.7 

0.32l^-f 

4.6 ± 1.7 

24.0 

176.4 ± 13.4 

199.8 ± 13.4 

70.0 

39.0 ± 28.3 

175.4 ± 28.3 


^ Color corrections for spectral flux density measurements were calculated from methods described in the IRAC and MIPS handbooks. 
^ The 5.8 /am coma photometry measurement is most likely contaminated by presence of the diffraction spike, leading to the higher 
measurement than the 5.8 nucleus photometry measurement. 


4. Dust Compositions 

4 . 1 . Extended Souree Calibration and Speetral Extraetion 

IRS spectra are calibrated with spectroscopic observations of point sources, as we have explained in several prior papers 
(Lisse et ah 2006, 2007; Sitko et ah 2011). The IRS slits are narrow with respect to the PSF of the telescope, i.e., they 
do not encompass 100% of the PSF at any wavelength and the fraction of the PSF encompassed varies with wavelength. 
Because the IRS instrument is calibrated with observations of spatially unresolved stars, the slit-losses may be ignored 
when working with calibrated observations of point sources. Ignoring the slit-loss does, however, affect the spectral shape 
of calibrated observations of sources larger than the PSF of the telescope (i.e., like what we have for SWl). We used the 
photometry from the IRAC 8.0 /im and MIPS 24.0 jam images, with an aperture of 9".0 radius, to tie the spectral flux 
density values with the IRS spectrum, which can be seen in Figure 13. 

4 . 2 . Speetral Modeling 

The infrared emission from a collection of dust near a comet nucleus has been explained in prior articles (Lisse et al. 
1998, 2004, 2005, 2006, 2007; Reach et al. 2010; Sitko et al. 2011) and this procedure of spectral analysis is used here. 
The thermal emission from the dust coma is given by 

^A.mod = ^2 

where T is the temperature for a particle of radius a and composition i at heliocentric distance r/j,, at a distance A from 
the observer, Bx is the blackbody radiance at wavelength A, is the emissivity (emission efficiency) of the particle 

of composition i at wavelength A, dnjda is the differential particle size distribution (PSD) of the emitted dust, and the 
sum is over all species of material and all sizes of particles for the dust. Spectral analysis consists of calculating the 
emission flux for a model collection of dust, and comparing the calculated flux to the observed flux. The emitted flux 
depends on the composition (location of spectral features), particle size (feature to continuum contrast), and the particle 
temperature (relative strength of short vs. long wavelength features). Prior to the modeling, the effects due to reflected 
sunlight and nucleus thermal emission need to be checked, and if necessary, removed. The IRS spectrum was divided by a 
Planck function, with temperature T = 140 K to remove the underlying blackbody profile of the dust grains. We excluded 
data below 8.0 /im due to the noise below this wavelength. The resulting relative emissivity spectrum along with a six 
component spectral model are shown in Figure 14. 


The model and its component-by-component decomposition is shown in Table 3 and Figure 15. Note that while we 
limit ourselves to six components, it is possible to lower the best-fit by adding components. It may also be possible to 
lower the best-fitting by choosing a different set of six components. For these reasons, we emphasize that our result 
here for the dust composition is not necessarily unique. However we claim that the structure of the spectrum indicated 
that the six components we present in Figure 15 have at least a reasonable likelihood of actually being in the grains. 
Specifically: The amorphous olivine is clearly seen in the spectrum by the broad emission feature from 15-20 /im. There 
is some evidence also for it in the 9-12 /im emission feature although the spectral data are of poor signal-to-noise ratio 
there that we do not want to rely on that. Water ice is revealed in the spectrums the increase in emissivity past 25 /im. 
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Figure 12: Shown are the photometry measurements from Spitzer observations. Curves shown are of the NEATM fit to 
the measurements and reflected solar flux. 1 -a error bars are included on the spectral flux density measurements. 



Figure 13: Shown are the coma photometry measurements for the IRAC and MIPS images for a 9".0 radius aperture, 
along with the Spitzer IRS spectrum. 
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Figure 14: IRS relative emissivity spectrum of SWl and best-fit six component spectral model. Blue points show the 
extracted data from Spitzer IRS observations, along with green error bars for each data point. The reduced = 1-036 
for model fit. The spectral decomposition of this model is shown in Figure 15. The model does a reasonable job of 
reproducing most of the broadband features, but there are wavelength regions (e.g. 12-16 /im) where the model deviates 
noticeably from the IRS observations. 


Table 3 


Composition of the Best-fit Model^ to the Spitzer IRS SWl Spectrum 


Species 

Weighted 
Surface Area 

Density 
(g cm“^) 

M.W. 

^moles 

(Relative) 

Model Tmax^ 
(K) 

Model 

(if not included) 

Amorph Olivine (MgFeSi 04 ) 

0.515 

3.6 

172 

1.08 

140 

62.10 

Forsterite (Mg 2 Si 04 ) 

0.210 

3.2 

140 

0.48 

140 

5.33 

Diopside (CaMgSi206) 

0.110 

3.3 

216 

0.17 

140 

2.65 

Pyrrhotite (MgFeS) 

0.140 

4.5 

84 

0.75 

140 

6.25 

Amoprh Carbon (C) 

0.675 

2.5 

12 

14.1 

190 

2.91 

Water ice (H 2 O) 

0.700 

1.0 

18 

3.89 

92 

36.61 


^ Best-fit 6-component model with power-law particle size distribution dn/da ^ a ^ (Note: using 40:60 Mg:Fe Amorph Olivine) . 

^ ^moles^ ^ density(i)/molecular weight(i) ^normalized surface area(i). Errors are ± 10% for all species, except Amorph Carbon, which has ± 20%. 
^ Temperature of submicron- to micron-sized dust grains in K. Errors in temperature of dust are ± 5 K (2-cr). 
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Figure 15: Best-fit six component spectral model is shown with the individual six component lab spectra used for fitting. 
The model was built on a linear combination of the six individual lab spectra, varying the linear expansion coefficients 
until arriving at a minimum deviation from the IRS relative emissivity. When comparing the individual spectra with the 
IRS emissivity in Figure 14, the signatures for the compositional elements become noticeable (e.g. the broad absorption 
feature of amorphous olivine around 13.5 /im or the narrow features of forsterite from 17-34 /rm). 


Diopside and forsterite are known by several emission features that seem to match such features in the spectrum itself 
long ward of 12 fim. Again, we do not want to rely on the features from those components at shorter wavelengths because 
the data are so uncertain. Pyrrhotite manifests itself as the broad emission around ^24-33 /im and it has the effect of 
bumping up the entire spectrum at those wavelengths. Finally, amorphous carbon helps the model fitting short ward of 
12 /im due to its upturn in flux relative to a 140 K blackbody there, but the high uncertainty of the data in this region 
makes it impossible to be definitive. 

The derived dust temperature of 140 +/- 5 (2-cr) K is close to what is expected for submicron- to micron-sized Mg-Fe 
amorphous silicate and amorphous carbon grains at ^ 5.2 to 5.6 AU from the Sun (e.g.. Banner et al. 1997). The derived 
particle size distribution (PSD) from the linear emissivity model is dnjda ^ ^- 3 . 75 ±o.i 5 (2-cr), indicating a comet with a 
high relative abundance of small dust particles. The dust is amorphous-dominated, which has been seen in only one other 
(bright) comet, 73P (Sitko et al. 2011 ). Amorphous, primitive dust is unusual, and likely indicates that the material 
composing the refractory portion of SWl was never heated to high temperatures, which is in contrast to crystalline silicate 
signatures observed in other comet dust for which we have good spectra (like 9P/Temple 1, Hale-Bopp (C/1995 01), etc.). 

Mie models of the same Spitzer spectrum shows similar results, with a somewhat steeper PSD (dn/da ^ q^- 3 . 9 ±o. 20 ^ 
but the difference between this and dn/da ^ q-3.75±o. 15 within the normal range of modeling error). These results 
strongly suggest that the most likely Mie models are those with a high relative abundance of amorphous carbon, leading 
to a high “superheat”, large T(dust)/T(blackbody) (Gehrz and Ney 1992), and low average dust albedo, as amorphous 
carbon grains are very efficient absorbers of optical solar energy (i.e., they are very dark), but only moderately good at 
re-radiating their energy in the thermal infrared. 


5. Conclusions and Discussion 

The reanalysis of archival Spitzer thermal observations of SWl has led to new measurements of physical properties 
of the nucleus. The effective radius measurement of R = 30.2 ^3 7 km is within the 1 -cr uncertainties of several recent 
measurements (Stansberry et al. 2004; Stansberry et al. 2007), signifying a possible honing of the nuclear size measure¬ 
ment. The Hv implied by the newly derived radius suggests a lower albedo for SWl. While the choice of geometric albedo 
and our determined value of the radius affects the By value, we would like to emphasize that our choice of a geometric 
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Figure 16: Histograms of measured beaming parameters from current journal articles (see conclusion for list of beaming 
parameter sources). The vertical bar in each plot signifies our measured value of 77 = 0.99 for SWl. The left plot shows the 
distribution of 57 JFC and SWl’s measured value is towards the mean of the distribution. All JFC r] values were acquired 
from Spitzer observations in SEPPCoN (Survey of the Ensemble Physical Properties of Comet ary Nuclei) Eernandez et 
al. (2013). The middle plot shows SWl’s placement in the ensemble of p values for 58 Centaurs. Again SWl’s measured 
value places it towards the middle of the distribution. Histogram colors represent the telescope used for observations 
(blue:Spitzer (Eernandez et al. 2013; Stansberry et al. 2008); green: WISE (Bauer et al. 2013); red:Herschel (Duffard et al. 
2014; Eornasier et al. 2013; Lellouch et al. 2013; Mommert et al. 2012 ; Pal et al. 2012 ; Santos-Sanz et al. 2012 ; Vilenius 
et al. 2012 )). Plot to the right is similar, but with SWl’s placement in the distribution of measured TNO p values. Again 
SWl’s measured value is towards the center of the distribution. 


albedo pv = 0.04 does not have much of an effect on the radius derived from thermal photometry. The following radius 
measurements were arrived at when varying the geometric albedo: (pv = 0.06, R = 30.2 km), (pv = 0.08, R = 30.2 km), 
and (pv = 0 . 2 , R = 30.1 km). All of the radius determinations are within the uncertainty of the value arrived at using pv 
= 0.04. If we choose to use pv = 0.04 then our measured radius leads to a value of Hy = 10.24. Stansberry et al. (2004) 
used a value of Hy = 11.0 (from Degewij et al. 1981; Meech et al. 1993) which with our radius determination yields a 
value of py = 0.025, i.e., a lower albedo and a value that is consistent with some TNOs (Stansberry et al. 2008). The 
size of SWl places it on the larger end of the Centaur size distribution. Considering the activity levels of SWl at such 
a large heliocentric distance (e.g. Paganini et al. 2013) signifies a large supply of volatile material. If the orbit ever gets 
perturbed sending SWl into the inner solar system it most likely will be an impressive sight. 

An infrared beaming parameter of 77 = 0.99 ^0 26 the middle of the 77 distribution for other Centaurs. A 

collection of IR beaming parameter measurements for a database of 57 Jupiter family comets (JEC), 58 Centaurs, and 
75 trans-Neptunian objects (TNO) (Bauer et al. 2013; Duffard et al. 2014; Eernandez et al. 2013; Eornasier et al. 2013; 
Lellouch et al. 2013; Mommert et al. 2012 ; Pal et al. 2012 ; Santos-Sanz et al. 2012 ; Stansberry et al. 2008; Vilenius et 
al. 2012 ) is shown in Eigure 16. The SWl measurement is highlighted by a vertical bar on the histograms. The measured 
value of 77 = 0.99 for SWl is towards the middle of each of the three small-body populations, as shown an in Eigure 16. 
A value of 77 = 0.99 for SWl could signify a surface with no radical topography (i.e. no high temperature regions causing 
infrared beaming). Alternatively, SWl could have a low thermal inertia, resulting in little or no night side emission. 

While our size measurement is within uncertainties of several previous measurements, the infrared beaming parameter 
of each analysis is considerably different. Stansberry et al. (2004) used a value of 77 = 0.62. If we use their measured 
spectral flux density values, a fixed 77 = 0.99 value, and our model we arrive at a best fit radius value of 28.6 km. 
Unfortunately, a reduced = 17.53 signifies a poor modeling fit. This fit can be seen in Eigure 17a. A similar analysis 
was performed on the Stansberry et al. (2007) spectral flux density measurements. A best-fit radius value of R = 26.8 


15 







































km was found with a smaller reduced = 4.75 value. This fit can be seen in Figure 17b. These values are shown to not 
undervalue previous measurements, but to emphasize the robustness in the measurements found in this analysis. Using 
archival Spitzer data products and newer more robust coma removal techniques, we have a size measurement that is in 
agreement with previous measurements, solidifying SWl’s place as a large Centaur. More importantly, our modeling has 
resulted in a beaming parameter value close to the average of a sampling of JFCs, Centaurs and TNOs. 




(b) 


Figure 17: Modeling fits to the Stansberry et al. (2004) and (2008) spectral flux density measurements. 


We were unable to find constraints for the 5.8 fim geometric albedo of SWl with the thermal observations alone (no 
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coinciding observations of SWl in the visible were available). Thermal modeling resulted in a value of 0.5 ± 0.5 which 
shows a high degree of uncertainty. While we do not claim an infrared geometric albedo measurement of 0.5, we would like 
to point out albedo measurements in this wavelength region are lacking from literature. Infrared albedos may be higher 
than the usual 0.04 values found for a geometric albedo in the visible. Although the uncertainty in our measurement 
is high as well, we would like to show that this magnitude of infrared albedo is not physically out of the question. A 
normalized reflectivity gradient S' = 14.94 ± 1.09 [% (1000 A)“^] (Duffard et al. 2014) and a V-band albedo of 0.04 would 
yield an albedo of 0.31 ± 0.10 in the near-infrared. Note that the normalized reflectivity gradient is a simplification and 
assumes a linear relation between albedo change and wavelength, which is truly not the case for mid-infrared albedos. It 
is mentioned here to just show the reader that the albedos of Centaurs tend to increase towards the infrared and that a 
value of 0.04 is not valid at these longer wavelengths. The high uncertainty in the derived IR albedo also means that SWl 
may have a much lower IR albedo that 0.5. 

Using a linear combination of silicate, amorphous carbon, water ice, and metal sulfide laboratory emissivities, shown to 
successfully fit and explain the 9P/Tempel 1, Hale-Bopp (C/1995 01), 17P/Holmes, and 73P/Schwassmann-Wachmann 
3 comet dust spectra (Lisse et al. 2007; Reach et al. 2010; Sitko et al. 2011), we find the fit for the Spitzer SWl 
spectrum shown in Figure 14. Containing 281 spectral measurements, this spectrum is much more highly constraining 
on the possible dust composition and size models. Only fits containing appreciable amounts of silicates, water ice, and 
amorphous carbon fit the data well. 

To summarize: 

• The measured value of R = 30.2 km places SWl on the larger end of the Centaur size distribution. This value 
has been shown to be consistent (within 1 -a) of several previous measurements. 

• An infrared beaming parameter value, 77 = 0.99, measured for SWl is towards the middle of distributions of measured 
7 ] values for JFCs, Centaurs, and TNOs. 

• The 5.8 /rm infrared albedo was measured to be 0.5 ± 0.5, but with such a large uncertainty in this value it is still 
unclear how the reflectivity truly behaves in this wavelength regime. It was shown, however, that using a normalized 
reflectivity gradient for SWl and a standard visible albedo for comets that a high (> 0.30) albedo can be achieved. 

• Finally, using Spitzer IRS observations a dust composition model was found to show SWFs dust is dominated by 
amorphous minerals, with possible signatures for more thermally processed minerals. 
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